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ABSTRACT 

We present new i^-band echelle spectra, obtained with the NIRSPEC spectrograph at Keck 
II, for the massive star cluster "B" in the nearby dwarf irregular galaxy NGC 1569. From 
4^ . spectral synthesis and equivalent width measurements we obtain abundances and abundance 

O^' patterns. We derive an Fe abundance of [Fe/H]=— 0.63 ± 0.08, a super-solar [a/Fe] abundance 

ratio of +0.31 ± 0.09, and an O abundance of [O/H]=-0.29 ± 0.07. We also measure a 
^ ■ low ^^C/^^C ^ 5±1 isotopic ratio. Using archival imaging from the Advanced Camera 

0^ ■ for Surveys on board HST, we construct a colour-magnitude diagram (CMD) for the cluster 

Ci . in which we identify about 60 red supergiant (RSG) stars, consistent with the strong RSG 

features seen in the i7-band spectrum. The mean effective temperature of these RSGs, derived 
from their observed colours and weighted by their estimated i^-band luminosities, is 3790 
^ . K, in excellent agreement with our spectroscopic estimate of Teff = 3800 ± 200 K. From 

■ the CMD we derive an age of 15-25 Myr, slightly older than previous estimates based on 

\ integrated broad-band colours. We derive a radial velocity of (vr) = — 78 ± 3 km/s and a 

l/*) . velocity dispersion of 9.6 ± 0.3 km/s. In combination with an estimate of the half-light radius 

O : of 0'.'20 ± 0'!05 from the HST data, this leads to a dynamical mass of (4.4 ± 1.1) x 10^ Mq. 

The dynamical mass agrees very well with the mass predicted by simple stellar population 
models for a cluster of this age and luminosity, assuming a normal stellar IMF. The cluster 
I core radius appears smaller at longer wavelengths, as has previously been found in other 

. extragalactic young star clusters. 
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1 INTRODUCTION 

Massive star clusters are potentially useful test particles for con- 
straining the evolutionary histories of their host galaxies. They can 
remain observable for the entire lifetime of a galaxy (as illustrated 
by the old globular clusters which surround every major galaxy), 
and they are bright enough to be studied in detail well beyond the 
Local Group. The internal velocity broadening of their spectra is 
typically only a few km/s, making detailed abundance analysis at 
high spectral resolution feasible. In a previous paper we have taken 
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a first step towards exploiting this potential by analysing H and K- 
band spectra of a young massive star cluste r (YMC) in the nearby 
spiral galaxy NGC 6946 jLarsen et al .120061) . Here we apply a sim- 
ilar analysis to one of the young star clusters in the nearby (post-) 
starburst galaxy NGC 1569. 

NGC 1569 was one of the first galaxies in which the pres- 
ence of exceptionally bright, young star clusters was suspected. 
A spectrum of one of the two brig ht "stellar con densations" in 
NGC 1569 was obtained already by iMavalll (Il935h . although the 
true nature of these objects was probably first discussed in de- 
tail by lArp & Sandagd (Il985l) . They found that the spectra are 
of composite nature, and also noted that observations with the 
Hubble Space Telescope (HST) would definitively settle the is- 
sue whether or not these objects are really "super star clus- 
ters". Indeed, pre-refurbishment mis sion HST observations by 
lO'ConnelL Gallagher & Hume3 (Il994h settled the issue by show- 
ing that both objects are extended, with half-light radii of about 
two pc. Based on high-dispersion spectroscopy from the HIRES 
spectrograph on the Keck I telescope, a dynamical mass of ^ 10^ 
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M0 was soon after estimated for the brighter of the two objects, 
NGC 1569-A (IHo & FiHpDenkolll996l : lsternberdl 19981 ). 

The clusters in NGC 1569 are prime candidates for abundance 
analysis in the near-infrared where they are very bright. The IR 
holds a particular advantage over optical studies for NGC 1569 due 
to the large amount of foreground extinction. [Gilbert & GrahamI 
(l2002h found the integrated if -band spectra of both clusters A 
and B to be well approximated by red supergiant templates. Al- 
though NGC 1569-A is the brighter of the two, it is actually a bi- 
nary cluster itself with some e vidence for a (sniall) age difference 
between the two c omponents fa e Marchi et al.l [19971 : iMaoz et al.l 
l200ll : IOriglia et al . 2001). Whether or not the two components are 
physically connected or the result of a chance projection is un- 
known. In this paper we concentrate on NGC 1569-B. 

In addition to providing insight into the histories of their host 
galaxies, young star clusters with masses in excess of 10^ Mo also 
offer an excellent opportunity to study large samples of coeval mas- 
sive stars. Such clusters are rare in the Milky Way and even in the 
Local Group, so it is necessary to extend the search to a larger vol- 
ume. The most massive kn own young star cl usters in the Milky 
Way disk are Westerlund 1 fclark et al." 2005) and an object dis- 
covered in the 2MASS survey ( Figer et al. 2006), both of which 
may have masses approaching 10^ M©. Young clu sters of similar 
masses are found in the Large Magellanic Cloud (Ivan den Berghl 
I1999I) . but these pale in comparison with objects like NGC 1569-A 
and NGC 1569-B. 

We have obtained new NIRSPEC if -band spectra of 
NGC 1569-B, optimised for abundance analysis, and we addition- 
ally present photometry for individual stars in the cluster derived 
from archival observations with the high resolution channel (HRC) 
of the Advanced Camera for Surveys (ACS) on HST. The HST data 
provide an independent verification of the stellar parameters (no- 
tably Teff and log g) derived from the spectral analysis, and also al- 
low us to compare the observed colour-magnitude diagram (CMD) 
with standard isochrones. Finally, we use new measurements of the 
structural parameters and velocity dispersion to derive the cluster 
mass and mass-to-light ratio and compare with predictions by sim- 
ple stellar population (SSP) models. 

We begin by briefly describing the data in 32 The analysis and 
our main results then follow in 321 where we first discuss the near- 
infrared spectroscopy and the abundance analysis (ED- We then 
proceed to construct a CMD from the HST imaging ( ^3.21 ) from 
which we derive stellar parameters for the red supergiants (RSGs) 
that are compared against the spectroscopic results ( ^3.2.31 ) and the- 
oretical isochrones ( ^3.2.41 ). In ^3.3l we measure structural param- 
eters for NGC 1569-B and combine these with velocity dispersion 
measurements to derive a dynamical mass estimate. Finally, some 
additional discussion and a summary are given in 33] 



2 DATA 

2.1 NIRSPEC data 

We obser ved NGC 1569-B in t he if -band with the NIRSPEC spec- 
trograph jMcLean et al.lll998b on the Keck II telescope on 2006 
November 8, using the high resolution echelle mode with a slit 
width of OV432 (3 pixels). This yields a spectral resolution of 
A/ A A = 25000. Since the NIRSPEC detector does not cover 
the entire ii-band echellogram in a single exposure, we optimised 
the tilts of the NIRSPEC gratings so that the number of features 
suitable for abundance determination was maximised. To facilitate 



easy sky subtraction, we followed the standard strategy of obtaining 
pairs of exposures nodded a few arcseconds along the slit. A total 
of 20 integrations, each with an exposure time of 300 s, were made, 
resulting in a total exposure time of 6000 s. One-dimensional spec- 
tra were extracted from each pair of nodded exposures and wave- 
length calibrated using the REDSPEC IDL package written by L. 
Prato, S. S. Kim & I. S. McLean. The individual 1-D spectra were 
then co-added. The final summed 1-D spectra have a S/N of about 
100 per pixel, or 175 per 3 pixel resolution element. 



2.2 HST/ACS data 

We use observations of NGC 1569 obtained with the ACS HRC as 
part of programme GTO-9300 (PL: H. Ford). The HRC has a field 
of view of about 29" x 2&' and a pixel scale of - 0V027 pixel" \ 
but with significant variations across the field due to the geometric 
distortions in ACS. For the photometry of individual stars we use 
the F555W {'^ V-band) and F814W (?^ I-band) images, each ex- 
posed for 3 X 130 s. For measurements of structural parameters we 
also include F330W U-band) data (3 x 220 s). The data were 
downloaded from the HST archive at the Space Telescope - Eu- 
ropean Coordinating Facility (ST-ECF) and ind ividual exposures 
were co-added with the MULTIDRIZZLE task ( Koekemoer et al.l 
l2002h . Multidrizzle filters out cosmic ray hits and detector defects, 
corrects for geometric distortions and projects the images on an 
orthogonal pixel grid with a scale of exactly 0V025 pixel" ^. At a 
distance of 2.2 Mpc f ^2.3l ), this corresponds to a linear scale of 0.27 
pc pixel" ^. 



2.3 Distance and Reddening of NGC 1569 

While it is clear from the high degree of resolution into individual 
stars that NGC 1569 must be nearby, the exact distance remains 
quite uncertain. One uncertainty is in the correction for foreground 
extinction, for which the NA SA Extragalactic Datab ase (NED) 
gives eithe r Ab ^ 2.03 mag iBurst ein & Heiles"l982') or Ab ^ 
3.02 mag ( Schlegel et aDll998h . Ilsrael (1988) estimated E{B - 
V) = 0.56±0.10orAs = 2.30d=0.40 mag (for =4.1). Com- 
bined w ith the apparent distanc e modulus (m — M)b — 29.0 de- 
rived bv lArp & Sandagd (Il985h . this leads to (m - M)o 26.7 =b 
0.6 or D — 2.2 zb 0.6 Mpc, which is the value used bv'Origlia et al.l 
(l200lh and several other studies in the literature. Based on photom- 
etry for a small number o f yellow and blue supergiants and also 
using E{B -V) ^ 0.56 JO'ConnelL Gallagher & Huntej 119941) 
estimated D — 2.5 =b 0.5 Mpc. Using observati ons of the red gi- 
ant branch tip,'Makarova & Karache ntsevI (12003 ^ derived two pos- 
sible distances of either D = 1.95 ± 0.2 Mpc or D = 2.8 ± 0.2 
Mpc, both assuming the Schlegel et al. foreground extinction. For 
the lower .Israel (1988) extinction, the distances would instead be 
D ^2.1 ± 0.3 Mpc or D = 3.9 =b 0.3 Mpc. In th is paper we use 
the same value of 2.2 Mpc as Origlia et al. (2001") and the extinc- 
tion value from Israel ( 1988), i.e. ^(B-y) = 0.56 ov Ab = 2.30 
mag {Av = 1.74 mag). The abundance analysis is not affected by 
either of these parameters. For the analysis of the colour-magnitude 
diagram and the cluster mass we will give further comments when 
needed. 
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Table 1. Adopted stellar atmosphere parameters and abundance estimates for cluster B in NGC 1569. 



Teff [K] 


logg e[Kms-i] 


[Fe/H] 


[O/Fe] 


[Ca/Fe] 


[Si/Fe] 


[Mg/Fe] 


[Ti/Fe] 


[a/Fe]^ 


[Al/Fe] 


[C/Fe] 


3800 


0.0 3 


-0.63 


0.34 


0.33 


0.33 


0.23 


0.33 


0.31 


0.23 


-0.27 


±200 


±0.5 


±0.08 


±0.10 


±0.09 


±0.11 


±0.13 


±0.10 


±0.09 


±0.11 


±0.10 



[a/Fe] is the average [< Ca, Si, Mg, Ti > /Fe] abundance ratio. 

NGC1569 - SSC B 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




' ' ' ' ' ' ' ' 



' ■ ' ' I I I I 



1.0 



0.8 - 



0.6 
0.4 
0.2 



1.615 1.616 1.617 1.618 1.619 1.620 

T 



1.675 1.676 



1.677 






1 1 










- ' ^11 




i: 


Fe ¥ 




Fe - 


OH 1 
Si 














wavelength (/xm) 



Figure 1. Near-IR spectra of the SSC B of NGC 1569. Observed spectrum: open circles. SoHd line: best-fitting RSG spectrum. Dotted lines: Model spectra 
with abundances scaled by a factor of 0.5 and 2.0. A few atomic and molecular features of interest are also marked. 



3 ANALYSIS AND RESULTS 
3.1 Chemical abundances 

In general, abundance analysis from integrated spectra requires full 
spectral synthesis techniques to properly account for line blending 
and population synthesis to define the dominant contribution to the 
stellar luminosity. The near IR stellar continuum of young stellar 
clusters and st arburst galaxies is alm ost entirely due to luminous 
RSGs (see e.g. lOriglia & Olivall200Ql for a review) and usually it 
dominates over nebular and dust emission. This represents a major, 



conceptual simplification in population and spectral synthesis tech- 
niques, making the interpretation of integrated spectra much eas- 
ier. The spectra can be modelled with an equivalent, average star, 
whose stellar parameters (temperature Tgff , gravity logg and micro- 
turbulence velocity mainly depend on the age and metallicity. 

Our analysis of the NIRSPEC spectra of NGC 15 69-B closely 
follow ed that in our study of a YMC in NGC 6946 jLarsen et all 
l2006h . A grid of synthetic RSG spectra for different input atmo- 
spheri c parameters and abundances were computed, u sing an up- 
dated dOriglia. Rich & Castroll2002l : IOriglia et alJl2003h version of 
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the code described inlOrig lia. Moorwood & Olival (Il993h . Briefly, 
the code uses the LTE approximation and is based on molecular 
blanketed model atmospheres of Johnson, Bernat & Krupp (1980|) 
at temp eratures <4000 K. Re cently, the NextGen model atmo- 
spheres jHauschildt et al.ll999h have been also implemented within 
the code and tested, yielding very similar results. The code also 
includes several thousands of near IR atomic lines and molec- 
ular roto-vibrational transitions due to CO, OH and CN. Three 
main compilations of atomic oscillator stren gths are used, namely 
the Ku rucz databas and those publish ed bv lBiemont & Grevessd 
(Il973h and Melendez & Barbuv ("19991). For further details we re- 
fer to [Origha, Moorwood & Oliva ( 1993|); lOriglia. Rich & Castrol 
(l2002h and Origha et al. ( 2003). 

As was done for the YMC in NGC 6946, we estimated the ran- 
dom errors on the derived quantities using a number of test models. 
With respect to the best-fitting spectrum, we varied the temperature 
by ATeff of ±200K, the microturbulence by of =Fl-0 km s"\ 
gravity by A log g of ±0.5 dex and abundances varying accord- 
ingly by ^ ±0.2 dex, in order to still reproduce the depth of the 
observed features. By comparing the mean difference {5) between 
the observed and best-fitting model spectrum with the correspond- 
ing ^- value for each test model, we found that the te st models are 
at best statistically significant at > 2a level, only. See lLarsen et al.l 
12006) for details. 

Fig. [T] shows several regions of the cluster spectrum with the 
most prominent features indicated. As noted in several previous 
studies, the spectrum shows many absorption features indicating 
the presence of RSGs. The open circles are the observed spectrum, 
the solid line the best fitting RSG model spectrum and the dotted 
lines are models with a factor of two different abundances from the 
best fit solution. Synthetic spectra with lower element abundances 
are systematically shallower than the best-fit solution, while the op- 
posite occurs when higher abundances are adopted. 

By fitting the full observed iif-band spectrum and by mea- 
suring the equivalent widths of selected lines, we obtained the 
stellar parameters and abundance pat terns listed in Table [j] Ref - 
erence Solar abundances are from iGrevesse & Sauvall (Il998h . 
Among the measured quantities are an Fe abundance of [Fe/H] 
= -0.63 ± 0.08 and [0/H] = -0.29 ± 0.07 (i.e. 12 + log 
(0/H) = 8.54), which combine to [0/Fe] = 0.34 ± 0.10. We 
also measure an average [a/Fe]=0.31 ±0.09 enhancement for 
the other a-elements and ^^C/^^C ^ 5 ± 1. The latter ratio is 
even lower than the value we found for the YMC in NGC 6946 
(^^c/i^c = 8 ± 2). Isotopic ^^C/^^C abundance ratios below 
10 are quite common in Galactic RSGs lEamb ert & Tomkin 1974; 
iTomkin. Luck & Lamber j 1 19761 : iGonzalez & WallersteinI bOOOh 
and have also bee n found in the cluster NGC 330 in the Small 
Magellanic Cloud (iGonzalez & WallersteinI 19991) . The fact that the 
observed ^^C/^^C values in RSGs tend to be lower than standard 
stellar models predict ( Schaller et al. 1992) may indicate that extra 
mixing mechanisms are present. 

To our knowledge, no previous determination of the Fe abun- 
dance in NGC 1569 exists. Our value is s imilar to that of cool 
supergiants in the Small Magellanic Cloud (lHilllll999 ). Our oxy- 
gen abundance for NGC 1569-B, on the other hand, is more 
similar to that of young stars in the Large Magellanic Cloud, 
although the number of stars with reli able abundance determi- 
nations is small CRolleston et "Zl l2002h . Our O abundance for 
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Figure 2. Left: ACS/HRC F814W image of NGC 1569-B. The circles in- 
dicate the region for which a CMD was constructed. Right: The residual 
image produced by DAOPHOT/ALLFRAME. 



NGC 1569-B is about a factor of two higher than the commonly 
cited value of 12 + log (0/H) = 8.19 ± 0.02 for Hll regions in 
NGC 1569 dKobulnickv & SkillmanI 19971). Other measurements of 
the O abundance in NGC 1569 come from'Storchi-Bergmann et al.] 
(1994) who derived 12 ± Ior (0/H) = 8.37, intermedi ate be- 
tween our va l ue and that of iKobulnickv & SkillmanI (Il997l) . while 
iDevost et al.l (Il997h got 12 + log (0/H) = 8.26. All these previ- 
ous measurements refer to the ionised gas, however, and we are 
not aware of any previous measurement of the O abundance in 
NGC 1569 from stellar spectroscopy. 

From the spectral fits we also measure a heliocentric radial 
velocity of (^;r) = — 78 ± 3 km/s and velocity dispersion (cor- 
rected for instrumental broadening) of ~ 9.9 ± 1 km/s. This 
velocity dispersion is hig her than the 7.5 ± 0.2 km/s found by 
[Gilbert & GrahamI (l2002h . but agrees well with our measurement 
based on cross-correlation analysis ( ^3.3.41 ). The radial velocity it- 
self is in excellent agreement with the optical radial velocity of 
NGC 1569 listed in the Third Reference Catalogue of Bright Galax- 
ies (-74 ± 17 km/s; de Va ucouleurs et al.lll99lh . 

3.2 The duster colour-magnitude diagram 

3.2.1 Photometry 

The HRC images show NGC 1569-B resolved into its brightest 
constituent stars. This presents us with a welcome opportunity to 
check the properties of RSGs derived from the spectral fits. To 
study the CMD, we therefore carried out point- spread-function 
(PSF)-fitting photometry with the DAOPHOT III photometry pack- 
age (St etson 1987, 1990) . Stars were first detected in the F814W 
image with the FIND task, aperture photometry was obtained for 
the F555W and F814W images with the PHOT task and PSFs were 
constructed using 7 isolated stars scattered across the HRC detec- 
tor. A first iteration of PSF-fitting photometry was then carried out 
by running the ALLFRAME task on the F555W and F814W frames 
simultaneously. The FIND routine was applied a second time to the 
star- subtracted F814W image produced by ALLFRAME, and im- 
proved PSFs were constructed by running the PSF task once again 
on images where all stars detected in the first pass (except the PSF 
stars) had been subtracted. The ALLFRAME PSF-fitting was then 
repeated with the combined coordinate list. We found that two such 
iterations were sufficient to produce very clean final star- subtracted 
images. Fig.|2]shows the F814W image of NGC 1569-B in the left- 
hand panel and the residual image produced by the second ALL- 
FRAME iteration in the right-hand panel. The stars subtract out 
nicely without any systematic residuals, even quite close to the cen- 
tre of the cluster, indicating the PSF-fitting procedure went well. 
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Figure 3. Colour-magnitude diagram for NGC 1569-B. The three curves 
show the 75% (dashed line, red in the on-line edition), 50% (solid line) and 
25% (dotted-dashed line, blue in the on-line edition) completeness limits 
based on artificial star experiments. 



Figure 4. Colour-magnitude diagram for artificial objects after recovery 
and PSF-fitting photometry. Objects with input magnitudes of mpsssw = 
20, 21, ...24 and mpsssw — ^F814W = 0? 1? 2, 3, 4 are shown here. 



For calibration to the VEGAMAG system we used the photo- 
metric zero-points for the HRC published on the ACS world- wide- 
web-page (2;f555W = 25.255 mag, 2;f814w = 24.849 mag). To 
fix the zero-points of the PSF fitting photometry we carried out 
aperture photometry in a radius of 5 pixels (0^.^125) for a number of 
bright stars on images where all other stars had been subtracted and 
calculated the offset between the aperture and PSF photometry for 
these stars. We estimate that the uncertainty in tying the PSF and 
aperture photometry together is about 0.01 - 0.02 mag. The aper- 
ture corrections from 5 pixels to infinity were take n from the tables 
of encircled energy distributions in lSirianni et al.l(l2005l) . 



3.2.2 Completeness and Contamination 

The resulting CMD is shown in Fig.[3]for the radial range 15 pixels 
< R < 50 pixels (0'.'37 < R < 1'.'25). Closer to the centre the 
crowding becomes too severe to obtain useful photometry and fur- 
ther out contamination from the general field population becomes 
an increa sing worry. Our CM D is consistent with the one already 
shown bv lHunter et al.! (l2000h , who used data from the Wide Field 
Planetary Camera 2, in containing a large number of RSGs. 

The severe crowding requires that careful attention be paid to 
completeness effects. We quantified the detection incompleteness 
and photometric errors by adding artificial stars in the same an- 
nulus for which the CMD was constructed. A common procedure 
when carrying out such artificial object tests is to add stars at ran- 
dom positions within an image; however, it is not necessary that the 
positions are truly random. As long as the underlying distribution of 
stars is random itself, no loss of generality results from adding the 
artificial stars in a semi-regular grid. The advantage is that they can 
be more closely packed this way. The only caveat is that care must 
be taken to sample different positions of the PSF at the sub-pixel 
level, in particular for undersampled images such as those from 
HST. We distributed the artificial stars in concentric rings with a 
radial separation of 10 pixels between each ring and a separation 
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Figure 5. Statistical decontamination of the CMD for NGC 1569-B. Top: 
the CMD for the annulus illustrated in Fig |2] Centre: the CMD for an an- 
nulus of larger radius, but identical surface area. Bottom: sources in the 
central panel which have a match within 0.5 mag in the top panel have been 
subtracted from the top panel. 
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of 10 pixels between stars within each ring. This approach auto- 
matically provides the required (pseudo-) randomness of the PSF 
positions at sub-pixel level and ensures that no two artificial stars 
are closer together than 10 pixels, thus avoiding "self-crowding" of 
the artificial stars. 

The artificial stars were recovered using exactly the same pro- 
cedure as for the regular photometry. In order to sample the incom- 
pleteness in different regions of the CMD, artificial stars with mag- 
nitudes between mpsssw = 20.0 and mpsssw = 26.0 (in steps 
of 0.5 mag) and colours between mpsssw — mpsiw = 0.0 and 
^F555W — ^F8i4w = 4.0 (stcps of 1 mag) were added. By repeat- 
ing the procedure 5 times, each with a different origin of the grid, 
we could sample a total of 346 artificial stars at every magnitude 
and colour within the annulus used for the photometry, allowing us 
to quantify the detection completeness with good statistical accu- 
racy. 

Fig. 131 shows the CMD for the recovered objects. For clar- 
ity, every second magnitude step has been omitted, and only one 
of the 5 tests at each colour/magnitude is shown. A comparison 
with Fig. [3] shows that the concentration of red supergiants at 
^F555w — 'niF8i4w ~ 2.5 and mF555w ~ 23 is real and not 
just an artefact of incompleteness towards fainter magnitudes. If 
any significant number of red objects with mpsssw ~ 24 were 
present, they would be clearly detected. On the other hand, in the 
blue part of the CMD the fainter magnitude limit appears to be pre- 
dominantly a detection limit. To put these statements on a firmer 
ground, we used the output from the completeness tests to calculate 
the magnitudes at which 25%, 50% and 75% of the artificial objects 
were recovered as a function of colour. These limits are shown in 
Fig. [3] and confirm that the concentration of RSGs is found well 
above any of these limits. The lower envelope of the remaining 
data points is consistent with the completeness limits inferred from 
the artificial star experiments. We conclude that the RSGs, which 
dominate the near-IR spectra, are well sampled by the ACS/HRC 
data, and we can use the photometry to constrain their properties 
further. 

The artificial star tests also show that much of the scatter in 
the colours of the RSGs can be accounted for by photometric er- 
rors. Because the errors on the mpsssw — ^f814w colours are 
dominated by the errors on mpsssw for red objects, the recov- 
ered artificial star magnitudes at the position of the RSGs scat- 
ter along a line which is nearly diagonal in the CMD. A similar 
trend is seen in the observed CMD in Fig. [3] This suggests that 
the RSGs in NGC 1569-B span a fairly narrow intrinsic range in 
^F555W — 'niF8i4w colour, and hence in Teff. 

Before proceeding, the question of contamination of the clus- 
ter CMD by the general field population of NGC 1569 must be 
addressed. Crowding prevents us from measuring stars closer than 
about 0^38 from the centre of the cluster, which is about twice the 
half-light radius. In the top panel of Fig. [5] we again show the CMD 
within 0^38 < R < 1^25, now corrected for Galactic extinction 
and converted to an absolute magnitude scale. The centre panel 
shows the CMD for a larger annulus with the same area. Complete- 
ness effects are expected to be less severe in this less crowded re- 
gion, yet the population of RSGs is virtually absent. This supports 
the notion that the RSGs in the top panel are i ndeed predominantl y 
associated with cluster B as also concluded bv lHunter et al.l (l2000l) . 
The bottom panel shows the statistical difference between the top 
and centre panels, constructed by removing every star in the central 
annulus which has a matching datapoint within 0.5 mag in m f555W 
and mF8i4w in the background annulus. For the RSGs there is no 
significant difference, although some fraction of the blue stars are 



removed. Since we have not applied separate completeness correc- 
tions in the two annuli, it is likely that too many blue stars have 
been subtracted near the detection limit. 

3.2.3 Comparison with spectroscopy 

Within the right-hand selection box in the bottom panel of Fig. \5\ 
we find 57 RSGs with a mean colour of (Mfsssw — Mf814w) — 
1.75 and (Mfsssw) = —5.6. For each star we can estimate the ef- 
fective temperature and bolometric corrections from the observed 
Mf555W — Mf814w colour, by using conversions computed with 
the code described by lOriglia & ^ Leithe rer (2000) and multiplying 
the model atmospheres bv lBessel et Z| imi) with the ACS filter 
responses. This yields (Tgff) = 3850 K, very close to the estimate 
of Teff = 3800 ± 200 K obtained from the spectral analysis. If 
we instead calculate the Teff values individually for each star and 
then compute an average weighted by the estimated if -band lu- 
minosity of each star, cooler stars will receive greater weights and 
the average Teff decreases to (Teff) = 3790 K, again in excellent 
agreement with the spectroscopic estimate. 

Different assumptions about the distance will not affect the 
photometric estimate of Teff, although different assumptions about 
extinc tion clearly will. If i nstead of Ab = 2.30 we assume Ab — 
3.02 (Schlegele taDll998h . then Teff increases by about 200 K. 

Having estimated the mean absolute magnitude and effective 
temperature of the RSGs from the CMD, the stellar diameter fol- 
lows and we can estimate the surface gravity if the mass is known. 
For an age of lo^ t — 7 .3 f ^3.2.4K isochrones from the Padua 
group (iGirardi et al.ll200Qh give a turn-off mass of 12.4 M©. As- 
suming Teff = 3800 K we get a bolometric correction of —1.4 
mag from the Kurucz modelo Then the mean RSG luminosity is 
LnsG /Lq — 4.9 X 10^ and the radius is i?RSG /Rq — 510, which 
leads to log 5' = 0.1 (in cgs units). This is again very close to the 
value log ^ = 0.0 derived from the RSG model spectral fit. 

3.2.4 Model comparison and the ratio of blue to red supergiants 

Models for massive stars are still fraught with significant uncer- 
tainties, partly because there is still only a limited amount of data 
available to compare with. CMDs for YMCs such as NGC 1569-B 
can help constrain the models. 

The mean effective temperature of RSGs predicted by models 
from the Padua and Geneva ( Lejeune & Schaerer 200.1 ) groups are 
listed in Table [2 for two different metallicities: Z — 0.004 and 
Z = 0.008. The Geneva models are available for "standard" and 
"enhanced" mass loss rates. Here we use the models with enhanced 
mass loss since these compare most favourably with the data. We 
have selected the RSGs in the same way as for the empirical data, 
i.e. RSGs are here defined as stars with Mfsssw — Mf814w > 
1.1. The (Teff) RSG values in Tableware the average of the Teff 
values listed in the isochrone tables for these stars, weighted by 
the luminosity in F 555W and the n umber of stars at each initial 
mass according to a iKroupal (l2002h IMF. Because all RSGs have 
very nearly the same mass for a given age, it makes virtually no 
difference what IMF is actually assumed and even for a flat IMF 
the mean effective temperature would change by only a few K. In 
all cases, the models predict the RSGs to be hotter than observed, 
with the Z = 0.008 Geneva isochrones coming closest to matching 
the observations in this respect. The Z — 0.008 Geneva models for 

^ Downloaded from |http ://kurucz.harvard.edu/grids .html] 
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Figure 6. Hess diagrams for Padua (left) and Geneva (right) Z = 0.008 isochrones for three different ages. The filled circles (red in the on-line edition) are 
the observed CMD for NGC 1569-B. 



standard mass loss produce RSGs which are about 100 K hotter 
than the enhanced mass loss models for the two youngest ages, 
while for the Z = 0.004 models (which are less affected by mass 
loss) there is no difference. 

In Fig. [6] and 17] we compare our background- subtracted CMD 
with the theoretical isochrones. Rather than simply drawing the 
isochrones as lines connecting the model points, we plot synthetic 
Hess diagrams which show the density of stars predicted by the 
isochrones. The CMDs confirm that the best match to the RSGs is 
obtained for Z = 0.008 with a slightly younger age being preferred 
for the Geneva models (logt ~ 7.2) than for the Padua models 
(lo g t 7.4). These ag es are older than the estimate of log t = 7.1 
by I Anders et al.l (|2004|) . based on integrated broad-band colours. 
However, if the distance of NGC 1569 is greater than we have 
assumed then the absolute magnitude scale will shift accordingly, 
favouring a younger age. It seems that the CMD favours a metal- 
licity for NGC 1569-B closer to Z = 0.008 than Z = 0.004, but 
we are reluctant to stress this point very strongly given the large 
uncertainties in the stellar models for these young ages. 

One specific "amplifier" of the uncertainties in models for 
massive stars is the number ratio of blue and red sup e rgiant s 
(BSG/RSG), as discussed in detail by danger & Maeded (Il995h . 
While it is relatively straight-forward to identify the RSGs in 
NGC 1569-B, it is, unfortunately, less clear how to select a sample 
of BSGs for comparison with the models as these are not clearly 
separated in the CMD. The Hess diagrams in Figs. [6] and |7] pro- 
vide some guidance and suggest that BSGs within the relevant age 



range should be brighter than My ~ —5.5. We therefore select 
BSG candidates within the left-hand box drawn in the bottom panel 
of Fig. [5] This results in 22 BSGs and 57 RSGs and thus a ratio of 
BSG/RSG = 0.39 ± 0.10. In Table [2 we list the predictions by 
the Padua and Geneva models, obtained by counting the BSGs and 
RSGs within the boxes shown in Fig. [6] and |7] We have used the 
same colour range as in Fig. [5j but adjusted the magnitude limits 
for each isochrone to avoid "contamination" from main sequence 
stars while making sure that the full range of BSG luminosities 
is covered in every case. Only the youngest Z = 0.008 Geneva 
isochrone can match the observed low BSG/RSG ratio, and in many 
of the other cases the discrepancy is very large, most notably for the 
metal-poor Geneva isochrones. 

More generous selection limits for BSGs in the observed 
CMD will reduce the discrepancy, but in any case it is clear that 
the isochrones do not reproduce the observed CMD in detail. The 
problem of the BSG/RSG ratio will be addressed in more detail in 
forthcoming pape rs, using a lar ger sample of star clusters of differ- 
ent ages (see also iLarsenlboOVh . 



3.3 Structural Parameters and Dynamical Mass 

3.3.1 Cluster size and its Wavelength Dependence 

Having discussed the stellar content of NGC 1569-B, we now turn 
to a reexamination of its structural parameters and dynamical mass. 
We measured the size of the cluster in the ACS/HRC images using 
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Figure 7. Same as Fig. [51 except that isochrones for Z = 0.004 are shown. 



our ISHAPE code (lLarsenlll999h . The cluster size is determined 
by convolving a series of analytic model profiles with the PSF and 
varying the mo del parameters until the Ijest fi t is obtained. We used 
"EPF" profiles (lElson. Fall & Freemanll 1987b of the form 

IiR) = Io[l + ir/rcf]~\ (1) 

where rc is the core radius and 77 determines the steepness of the 
profile at large radii. Because neither the cluster itself nor the fitting 
profile has a clearly defined outer boundary, we experimented with 
various fitting radii up to 50 pixels (15 pc). Fig. [8] shows the im- 
age of the cluster, the best-fitting model, and the residuals for each 
band, for the largest fitting radius of 50 pixels. The resolution into 
individual stars means that the xl values of the fits never get close 
to unity, but otherwise the models provide a good fit to the data. 

The PSFs used for the size measurements were based on 
DAOPHOT PSFs in the inner regions (R < 5 pixels). We know 
from the smooth residual images produced by the ALLFRAME 
photometry that these PSFs accurately represent the PSFs of our 
images. However, for the ISHAPE analysis we need to know the 
PSF at larger radii than for the photometry, and there the construc- 
tion of accurate empirical PSFs becomes limited by decreasing S/N 
and non-uniformity of the background. For the analysis of struc- 
tural parameters we therefore made use of composite PSFs, con- 
structed from the DAOPHOT PSFs for R < 5 pixels and using 
models computed by the TINYTIM package (iKrist & Hookll 19971) 
further out. 

Table[3]lists the Full Width at Half Maximum (FWHM) along 
the major axis, the minor/major axis ratio (B/A), the slope param- 



eter ?7, and the half-light radius Rh for various fitting radii. All fits 
consistently yield an axis ratio of B/A 0.8, confirming the conclu- 
sion from earlier work that NGC 1569 -B is significantly elongated 
dO'ConnelL Gallagher & Huntei]ll994l) . For ry < 1 the total lumi- 
nosity of the model profile is infinite, and Rh therefore undefined. 
For some fits 77 is close to unity, leading to unreasonably large Rh 
values. We list two values for Rh'. One assuming that the profile ex- 
tends to infinity (Rh{oo)) and another assuming that it is truncated 
at = 100 pixels (Rh{100)). For ry > 1 the two Rh values are 
similar, but for 77 1 they differ substantially. The cut at = 100 
pixels (30 pc) is chosen fairly arbitrarily, apart from the fact that 
it is very difficult to trace the cluster profile to greater radii in the 
images. 

Regardless of the fitting radius, all fits return a smaller FWHM 
at longer wavelengths. The formal errors on the FWHM values re- 
turned by ISHAPE are generally less than ±0.1 pixels, so the trend 
with wavelength appears significant. The FWHM is roughly equal 
to twice the core radius rc of the EFF profiles (for 77 = 1.0 this 
holds exactly) so the trend in FWHM with wavelength is basically 
a tren d in core radius. Suc h a trend seems consistent with the find- 
ing bv Hunter et al that NGC 1569-B has redder mpssew — 
^F555W colours near the centre. However, they also found that 
the mF555W — colours become bluer near the centre. 
In the top panel of Fig. |9] we show the mF33ow — mpsssw and 
^F555W — TnF8i4w colours vcrsus aperture radius for NGC 1569- 
B measured on the HRC data. The trends seen in this plot are very 
similar to those observed by Hunter et al., with mpssow — mpsssw 
becoming redder inwards and mpsssw — ^F814w getting bluer. 
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Table 2. Observed and model supergiant properties. BSG/RSG is the ratio 
of blue to red supergiants and (Tgff }rsg is the mean effective tempera- 
ture of red supergiants, weighted by their absolute MF555W magnitude. 
Red supergiants are here defined as stars having MF555W < —4.5 and 

^F555W - ^F814W > 1-1. 



Table 3. Structural parameters for NGC 1569-B. The FWHM and Rh are 
given in pixels (1 pixel = 0^025 = 0.27 pc (D / 2.2 Mpc)). 
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Figure 8. HRC images, ISHAPE model fits and residuals of NGC 1569-B 
in three passbands. These panels show a 100 x 100 pixels section, corre- 
sponding to 2'/ 5 X 2'/ 5. 



How can this be consistent with a monotonic decrease in cluster 
core radius with wavelength? 

The answer appears to be that a naive comparison of colour 
gradients, as in Fig. O neglects the effect of the PSF. Colour gra- 
dients can result from a combination of changes with wavelength 
in the structure of both the cluster and the PSF. The PSF behaviour 
with wavelength is not monotonic. Due to the "red halo" effect, 
more light is scattered to large radii at lon ger wavelengths, bu t 
the encircled energy distributions tabulated in lSirianni et alj (l2005h 
also show that more light is scattered to intermediate radii (0^.^ 1 - 



Table 4. A summary of properties of NGC 1569-B. is the half-light 
radius, mpsssw is our estimated total apparent magnitude in the F555W 
band, ar is the projected velocity dispersion, Vr the heliocentric radial ve- 
locity, Mvir the dynamical mass and ph the mean density within the half- 
mass radius. For quantities that depend on the assumed distance this depen- 
dency is specified explicitly. 



Rh (arcsec) 

^F555V^ (VEGAMAG) 
ar (km/s) 
Vr (km/s) 
Rh (pc) 

Mvir (Mq) 
Ph (Mq pC-3) 



0'.'20 ± 0'/05 
15.60 ±0.20 
9.6 ±0.3 
-78 ±3 
(2.1 ±0.5) (D/2.2Mpc) 
(4.4 ±1.1) 10^ (D/2.2Mpc) 
(2.5 ± 1.3) 10^ (D/2.2 Mpc)-^ 



1") in F330W relative to F555W. This is also illustrated by the fact 
that the aperture correction s to infinity display a minimum around 
5000A- 6000A (Fig. 10 in lSirianni et al.ll2005h . 

To investigate this further, we convolved the PSF in F330W, 
F555W and F814W with EFF models for the respective filters 
having FWHM of 5.0, 4.5 and 4.0 pixels and envelope slopes of 
?7 = 1.26, 1.23 and 1.20 (guided by the numbers in Table [3]). The 
centre panel of Fig. [9] shows the integrated colours versus aper- 
ture radius for these model profiles, normalised to match the ob- 
served colours at 0V5. Although the intrinsic cluster core radius 
decreases monotonically with wavelength by design, we reproduce 
the observed trend towards bluer mrsssw — mrsiw colours and 
redder mpssow — ^f555W colours near the centre. The bottom 
panel shows the colour- aperture relations for the pure EFF mod- 
els, as they would be observed in the absence of any convolu- 
tion with the PSF. In this case both the mpsssw — TnF8i4w and 
^F33ow — ^F555W colours bccomc redder near the centre. We 
conclude that the colour gradients observed bv lHunter et alJ (l200Qh 
(and also seen in the HRC data) are consistent with the monotonic 
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Figure 9. Top: mF330W - '^F555W and mpsssw - ^F814W colour 
versus aperture size for NGC 1569-B. To fit the two plots on the same scale, 
1 magnitude has been subtracted from mpsssw — ^F814W • Centre: Same, 
for the PSF in each filter convolved with EFF profiles of FWHM = 5.0, 4.5 
and 4.0 pixels in F330W, F555W and F814W. The model profiles have been 
normalised to the same colour as the observed profiles at a radius of 0(^5. 
Bottom: Same, but for the EFF profiles without any PSF convolution. 



decrease in the cluster FWHM with wavelength obtained from the 
ISHAPE fits. 

Because the envelope slope rj tends to be shallower at longer 
wavelengths, the correlation between FWHM and wavelength does 
not translate into a similarly strong trend of half-light radius Rh 
versus wavelength. The Rh{oo) values show a very larger scatter 
while most i^^(lOO) values are in the range 6-8 pixels or 0^.^15- 
0V20. IO'ConnelL Gallagher & Huntej (Il994h measur ed a half-light 
radius of Rh = 0^25 on HSTAVFPC images, while iHunter et all 
derived an almost identical value of Rh = 0^26 from 
HSTAVFPC2 images. Our measurement on the ACS/HRC im- 
ages is slightly smaller, but an increase in the cut-off radius by 
a factor of 2-3 would bring our Ru estimate into close agree- 
ment with previous studies. In the following, we therefore adopt 
Rh = 0V20 =b 0V05. 



3. 3. 2 Mass segregation ? 

NGC 1569-B is not the only YMC in which the s tructu re has been 
found to be wavelength dependent. iLarsen et"aD (l200lh noted that 
the colours of a YMC in NGC 6946 become bluer in larger aper- 
tures. A trend of smaller cluster size at longer wavelengths was also 
found for the cluster M82-F by McCra dv et all (l2005h . who argued 
that it might be du e to mass seg regation (although this result has 
been questioned bv lBastian et al.l 12007) after a detailed analysis of 
the extinction near M82-F). Because the RSGs are the most massive 
stars currently in the cluster, mass segregation would cause them to 
be more centrally concentrated, making the cluster appear smaller 
at longer wavelengths where the RSGs dominate. 

Mass segregation may be either primordial in nature (the most 
massive stars forming preferentially near the centre) or a result of 
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Figure 10. Integrated magnitude versus aperture radius for NGC 1569-B. 
The background was measured in an annulus between 2^0 and 3^0. 



dynamical evolution; in the latter case it is expected to occur on the 
half-mass relaxation time-scale which is 



1.7 X lQ^[i^/,(pc)]^/^Ari/2 



[m/Md 



11/2 



years 



(2) 



where N is the number of stars in the cluster and m is the mean 
mass per star (Spitzer 1987). For a total cluster mass of ~ 4 x 10^ 
Mq f ^3.3.5K a mean stellar mass of 0.6 M© (for a Kroupa IMF 
extending from 0.1 Mq to 100 M©) and Rh ^ 2 pc, Eq. ^ yields 
trh ~ 5x10^ years, about a factor of 25 longer than the current age 
of NGC 1569-B. Thus, the cluster is not expected to be fully mass 
segregated, but this statement is subject to a few caveats. Near the 
centre of the cluster, where the density is higher, the relaxation time 
will be shorter. Furthermore, massive stars can "sink" to the centre 
on shorter time scales due to energy equipartition. The timescale 
for this to occur is inversely proportional to the mass of the stars 
in question (Gerhard 2000), so stars with masses characteristic of 
the current RSGs might become segregated on ti me scales similar 
to the current age of NGC 1569-B. According to [McMillan et all 
(l2007h . dynamical mass segregation may occur even more rapidly if 
a cluster forms from subclumps which already have some degree of 
primordial mass segregation. Alternatively, mass segregation may 
be imprinted in the structure of the cluster from the beginning (e.g. 
references in lMcMillan et al.ll2007h . 

Whether or not mass segregation can actually explain the ob- 
served relatio ns between c luster structure and wavelength is an- 
other question. lFleck et aP 42006) calculated the difference Av-k 
between the y — K colour inside and outside the half-light radius 
for a mass-segregated cluster as a function of time and found this 
to be 0.05 - 0.15 mag. Using simple stellar population models by 
bruzual & Chariot! (l2003h we estimate Av-k/Av-i ~ 2, so we 
may expect mass segregation to affect Av-i by a few times 0.01 
mag. From the data plotted in the bottom panel of Fig.|9]we mea- 
sure a difference Av-i = 0.03 mag, which appears roughly con- 
sistent with the calculations by Fleck et al. So it seems that mass 
segregation might be responsible for the observed colour gradients, 
although it would be desirable to investigate this in more detail. 



3.3.3 Integrated magnitude 

For many of the same reasons that complicate measurements of the 
core- and half-light radius, the integrated magnitude of NGC 1569- 
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B remains uncertain an d a ran ge of numbers can be found in the 
literature. lAnders et all (l2004h list Mf555w — —12.83 for an 
R = 0^.^46 aperture. For their assumed distance of 2.2 Mpc, com- 
bined with reddening from Schlegel et al. and a correction for 
internal reddening in NGC 1569 of E{ B - V) = 005 mag, 
this corresponds to mpsssw = 16.30. lOriglia et"aD (l200lh list 
a reddening corrected apparent magnitude of mF555W,o = 14.5 
or mF555W = 16.24, close to the valu e given by A n ders et al. 
Using a larger aperture of = l'.'34, iHunter et al"] (|200Q) list 
Mf555W = -13.02 for E{B - V) ^ 0.56 and a distance of 
D — 2.5 Mpc, i.e. mpsssw — 15.70. 

Most of the differences between published photometry can be 
ascribed to different choices of aperture radii. Using the same aper- 
ture radii, we can largely reproduce the numb ers given above. For 
the apertures and background annuli used by I Anders et al.l (|2004|) 
and lHunter etaP (120001 we obtain mF555w(0'.'46) = 16.20 mag 
and mF555w( 1^.^34) = 15.66 mag. Figure [TO] shows the inte- 
grated apparent ttifsssw magnitude as a function of aperture ra- 
dius, demonstrating that the curve-of-growth never levels off com- 
pletely out to at least 2" . This is a generic difficulty for many young 
star clusters, which tend to have very extended outer envelopes 
with ?7-values close to 1 and hence poorly defined total luminosities 
(lLarsenll2004b . 

Per definition, the total luminosity of the cluster should be 
twice that contained within the effective radius. We can there- 
fore estimate the total integrated magnitude by measuring the 
magnitude within Rh and doubling the flux. From the data in 
Fig-Cniwe get mF555W = 16.79/16.55/16.40 mag in apertures of 
0'.'l5/0'.'20/0'.'25. Doubling the fluxes we then get mF555w(oo) ^ 
15.80 lb 0.20, where the quoted error is due to the uncertainty 
on the half-ligh t radius. This is consistent with the measurement 
of iHunteretal] (EoOO) but brighter than measurements made in 
smaller apertures. However, this estimate of the integrated mag- 
nitude may still be a lower limit, because we have used a half-light 
radius that was already corrected for broadening by the PSF. In re- 
ality, some fraction of the cluster light is scattered from the centre 
to larger radii by the PSF, and we therefore expect that an aper- 
ture with the same radius as the intrinsic Rh will contain slightly 
less than half of the total cluster luminosity. By carrying out aper- 
ture photometry on a simulated PSF, we estimate this effect to 
be about 0.20 mag, consistent with the encircled energy tables in 
ISirianni et al.l (l2005 h. Our final estimate of the integrated cluster 
magnitude is then mF555w(oo) ^ 15.60 ± 0.20. 

3.3.4 Velocity dispersion 

As noted in ^3.11 the best fit of our RSG model spectrum to the 
NIRSPEC data requires a further broadening of 9.9 ± 1 km/s af- 
ter accounting for the instrumental resolution. We have also de- 
rived the v elocity dispersion using the standard cross-co rrelation 
technique (iTonrv & Davislll979l : Iho & FilipDenkolll996h applied 
in many previous studi es. We refe r to our previous papers on this 
subject fordetails (e.g. L arsen et al. 2001; Larsen & Richtler 2004; 
[Uarsen et al.ll2004i) . We used two template stars, HR 2184 (spec- 
tral class Kill) and HR 2235 (F5 II). Fig. [TT] shows the cross- 
correlation function (CCF) for NGC 1569-B versus HR 2184 (thick 
curve) together with the CCF for the spectrum of HR 2235 (broad- 
ened by various amounts) versus that of HR 2184. The best match 
is obtained for a broadening of 10 km/s. 

In Fig. [12] we show the distribution of velocity broadenings 
derived from the cross-correlation analysis in the various orders of 
the NIRSPEC spectrum and the different combinations of template 
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Figure 11. Cross-correlation functions for order 48 of the NIRSPEC data. 
The CCF for NGC 1569-B versus HR2184 is shown as a thick line (blue 
in the on-line edition). Thin dashed lines show the CCF for HR2235 versus 
HR2184, with the spectrum of HR2235 being broadened by 1 km/s, 5 km/s, 
10 km/s, 15 km/s and 20 km/s. The narrowest CCF peaks correspond to the 
smallest amount of broadening. All CCFs have been shifted to a mean of 
km/s for easier comparison. 



stars. The best fits (highest CCF peaks) were obtained for orders 
45, 46 and 49 and the results for these orders are shown as a filled 
histogram. The mean is (Jr — 9.6 ±0.1 km/s. If we include orders 
47 and 48 the mean decreases to (Jr = 9.2 zb 0.1 km/s, still sig- 
nificantly greater than the 7.5 km/s measured bv lGilbert & GrahamI 
(2002). For the following discussion we will adopt Gr — 9.6 ± 0.3 
km/s, including the estimate of Or from the RSG model fit within 
our Icr error. 



3.3.5 Dynamical mass 

Assuming virial equilibrium and that mass traces light, the clus- 
ter mass is related to the 1-D velocity dispersion ar and half-light 
radius by 



(JrRh 



G 



(3) 



with a ^ 9.75. Inserting the Rh and Gr values derived in ^3.3.11 
and ^3.3.41 (summarised in Table |4]), we obtain a dynamical mass 
of Mvir = (4.4 zb 1.1) X lO^M© which is ab out a factor of two 
higher than the value of 2.3 x 10^ M© derived bv lGilbert & GrahamI 
(l2002h . The difference is mostly due to our greater velocity disper- 
sion. Note that Mvir scales linearly with Rh, and thus with the 
assumed distance D. Rather than trying to include an uncertainty 
on the distance, we have specified this dependency explicitly for the 
cluster mass listed in Table|4] This dynamical mass of NGC 1569-B 
is much greater than for any young star cluster currently known in 
the Milky Way and also well above the median mass of old globular 
clusters (though still comfortably within the mass range of Milky 
Way GCs). 

We can compare the dynamical mass estimate with the mass 
expected for a cluster of this age and luminosity, given a nor- 
mal stellar IMF. For Z = 0.008 and a Chabrier IMF, the 
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Figure 12. Distribution of velocity dispersion measurements from the 
cross-correlation analysis. The filled histogram shows the measurements 
based on the best orders (45, 46, 49) while the outlined histogram also in- 
cludes orders 47 and 48. The best orders yield a mean of (Jr = 9.6 ± 0.1 
km/s. 



iBmzual & Chariot! (l2003h SSP models predict a F-band absolute 
magnitude of My — 1.53 =b 0.15 mag per M© for logt = 
7.3 ± 0.1. For Mv ~ mrsssw - 26.7 - Ay ^ -12.85 ± 0.2 
the "photometric" mass is then (5.5 ±1.3) x 10^ M©, which is in 
agreement with the dynamical mass within the 1 a errors. 

The close agreement between the dynamical and photomet- 
ric mass estimates is probably somewhat fortuitous, given the un- 
certainties on the distance, reddening and age. In the absence of 
reddening, the dynamical-to-photometric mass ratio will scale as 
1/D. However, since the problem of determining the reddening 
and distance to NGC 1569-B are closely tied, the comparison is 
not as straight forward. To first order, the absolute magnitude of 
NGC 1569-B may be roughly independent of the reddening cor- 
rection, since most distance determinations to NGC 1569 are based 
on standard candles. The dynamical-to-photometric mass ratio will 
then scale as D instead of 1/D. All this neglects any possible de- 
pendency of the age estimate on distance and reddening. This will 
depend on the method used for age determination - ages based on 
integrated colours (or spectroscopy) will be distance independent 
and probably only weakly dependent on the reddening (since this 
is usually included as a free parameter in the fit) while CMD-based 
ages will depend on both. 

Another source of uncertainty is the parameter a in Eq. ([3]). 
In a mass segregated cluster, a will have a different value than we 
have assumed. Because of energy equipartition, the massive stars 
will acquire lower velocities and sink to the centre. Since measure- 
ments of the velocity dispersion and cluster size based on integrated 
light are dominated by these stars, t he ma ss derived from Eq. ^ 
will be underestimated. iFleck et al ](l2006^ model the evolution of 
a for dynamically mass segregated star clusters and give an analyt- 
ical approximation to the relative error Aa/a as a function of age 
and relaxation time. For an age of 20 Myr and trh — 500 Myrs, 
their equation (31) yields Aa/a ^ 13%. With this correction, the 
dynamical mass becomes Mvir,corr = (4.9 ±1.2) x 10^ M© for 



D — 2.2 Mpc, i.e. even in even closer agreement with the photo- 
metric estimate. 



4 DISCUSSION AND SUMMARY 

We have presented a detailed investigation of the stellar content 
and other properties of the massive star cluster 'B' in NGC 1569. 
Using new high S/N if -band echelle spectroscopy from the NIR- 
SPEC spectrograph on Keck II, we have carried out abundance 
analysis of red supergiant stars in the cluster. We find an iron 
abundance of [Fe/H] = -0.63 ± 0.08, close to that of SMC 
field stars. Our estimate of the oxygen abundance, [0/H] = 
—0.29 ± 0.07, is about a factor of two higher than that derived 
for Hll regions Istorchi-B er^mann et al ] |1994 iDevost et al.|[T997l : 
iKobuln ickv & Skillman 1997), and the resulting super-solar [0/Fe] 
abundance is unlike that observed in the Magellanic Clouds and 
young stellar populations in the Milky Way. However, according 
to our measurements NGC 1569-B also differs from these galaxies 
by having a higher [a/Fe] ratio. The difference between our [0/H] 
measurement for NGC 1569-B and the O abundance derived for 
Hll regions is greater than the formal uncertainties on either value, 
and our oxygen abundance is higher than the neb ular abundances 
in most dwarf galaxies (e.g. lVaduvescu et ani2007h . This raises the 
question whether one (or both) methods suffers from possible sys- 
tematic errors, or if there might be a genuine difference between 
cluster and Hll region abundances in NGC 1569 - perhaps due to 
an ab initio oxygen enhancement of the gas out of which the cluster 
formed. Since oxygen is the only element in common between the 
cluster and Hll region data, a detailed comparison of the abundance 
patterns is unfortunately not possible. However, it is of interest to 
note that the very strong Wolf-Rayet fe atures in the spec trum of 
cluster A also suggest a high metallicity (iMaoz et aDl200ll) . 

NGC 1569 is a popular target for testing models of chemi- 
cal evolution, and it is generally recognised that the strong galac- 
tic wind ob s erved in the galaxy must play an important role. 
iMartin et al ] (l2002h observed a super-solar [a/Fe] ratio for the 
wind, as we do for NGC 1569-B. However, it is unclear to what 
extent the two are related as material ejected in the outflow may 
not participate in star formation, and in any case not before it has 
had time to cool down and fall back into the galaxy. The chemi- 
cal evol ution of an NGC 15 69-type galaxy has been modelled in 
detail bv lRecchi etan ( l2006h . assuming a variety of bursty star for- 
mation histories. Their chemo-dynamical simulations show that el- 
ements produced in the last burst of star formation generally do 
not get mixed with the ISM in the galaxy, but instead get injected 
into the hot gas phase. None of their models produce an O abun- 
dance as high as the one observed by us, and the model that comes 
closest (12 ± log (0/H) ^ 8.4) severel y underpredicts the N/O 
abund ance of the Hll regions observed bv IKobulnickv & SkillmanI 
(Il997h . Generally, the best fits to the Hll region abundances are 
obtained for a "gasping" star formation history, but these models 
all predict a maximum 12 ± log (0/H) 8.1. Similarly, the O 
abundance predicted by the models of iRomano et al.l (l200d) for 
NGC 1569 never exceeds 12 ± log (0/H) = 8.41. We speculate 
that the galactic wind which is responsible for removing metals 
may be less efficient in doing so near the bottom of the potential 
well where the massive clusters formed. Unfortunately, none of the 
models include predictions for the wide range of other elements we 
are observing here. 

The Small Magellanic Cloud presents an interesting compar- 
ison case. Early spectroscopic studies and Stromgren photome- 
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try indicated that the young cluster NGC 33 is about 0.5 dex 
more metal-poor than the surrounding field jGrebel & Richtlej 
'l992, and references therein). High-dispersion spectroscopy by 
Hill (1999) showed a smaller and only marginally significant dif- 
ference between the field stars and NGC 330, although still in the 
sense that the cluster is more metal-poor than the field ([Fe/H] 
= -0.82 zb 0.10 vs. [Fe/H] = -0.69 dz 0.11). The [0/Fe] abun- 
dance ratios in these stars were all found to be sub- solar ([0/Fe] 
^ —0.15 to —0.3 dex), as in th e LMC and in young Galactic su- 
pergiants (iHill et al.lll997h . lHni 119991) also found the a-elements 
abundances (Mg, Ca, Ti) relative t o Fe t o be around Solar for 
NGC 330. iGonzalez & WallersteinI (Il999h found sUghtly lower 
iron abundanc es of [Fe/H ] = —0.94 zb 0.02 for K supergiants in 
NGC 330 thanlffii (Il999h and an [0/Fe] ratio closer to Solar. They 
derived somewhat enhanced [Ca/Fe], [Si/Fe] and [Mg/Fe] ratios 
(+0.18, +0.32 and +0.11), in contrast to sub-solar ratios derived 
for B stars. 

While there is some evidence for a difference in metallicity 
between NGC 330 and the SMC field, the situation there seems 
to contrast with the case of NGC 1569-B where we find a higher 
oxygen abundance of the cluster compared to the Hll regions and 
an enhanced [0/Fe] ratio. Our mean [a/Fe] = +0.31 ± 0.09 is also 
higher than that derived for the stars in NGC 330. One significant 
difference may be the mass of NGC 1569-B, which is probably an 
order of magnitude greater than that of NGC 330 (iFeast & Blackl 
ll98Qh . 

The spectral analysis returns a mean Teff = 3800 ± 200 K 
and log g ^ 0.0 for the RSGs in NGC 1569-B. We have checked 
these values using resolved photometry of the RSGs in the clus- 
ter, derived from archival HST/ACS images. About 60 RSGs are 
easily identifiable in the colour-magnitude diagram, and from their 
^F555W — ^F8i4w colours wc dcrivc a mean effective tempera- 
ture of (Teff) = 3850 K and log ^ = 0.1. An even closer match to 
the spectroscopic Teff is obtained if we weigh the Teff estimates for 
the individual stars by their if -band luminosities, in which case we 
get (Teff) = 3790 K. 

We have compared the CMD with Z = 0.004 and Z = 0.008 
isochrones from the Padua and Geneva groups. Both sets of models 
provide the best match to the observed colours and magnitudes of 
the RSGs for Z — 0.008 and ages of 15-25 Myrs, but no isochrone 
can reproduce the observed CMD in detail. For these sub-solar 
metallicities, the models predict higher (by up to a few 100 K) 
mean effective temperatures (Teff) for the RSGs than observed. 
Since the RSGs generally become cooler with increasing metallic- 
ity, models of higher metallicity are favoured, but models of Solar 
metallicity would be required to match the observed colours. For 
such models, however, the blue supergiants become much too cool 
to match the observations. The observed ratio of blue to red super- 
giants (BSG/RSG = 0.39 ± 0.10) is also significantly lower than 
most of the model predictions. Although we are detecting an im- 
pressive number of RSGs, it should be noted that since we are only 
resolving stars located well outside the half-light radius, the actual 
number of RSGs in NGC 1569-B is likely to be more than a factor 
of two greater than the number we have detected. 

We derive a velocity dispersion of 9.6 ± 0.3 km/s from the 
integrated spectrum o f NGC 1569-B, somewha t higher than the 
value of 7.5 km/s of "Gilbert & Graham| (l2002h . Combining this 
with an estimate of the cluster half-light radius of 0^20 ± 0V05 
measured on the ACS images, we obtain a dynamical mass of 
(4.4 ± 1.1) X 10^ ™^ excellent agreement 

with the mass predicted by simple stellar population models for a 
cluster of this age and luminosity and a standard IMF. A correction 



for mass segregation would bring the two estimates into even closer 
agreement. 

Our analysis of structural parameters reveals a decrease in 
core radius with wavelength for NGC 1569- B. This trend i s con- 
sistent with the colour gradients observed bv lHunteretaP feOOOl) 
and by us in the HRC data when variations in the PSF are taken 
into account. A similar colour gra dient was observed in an YMC 
in NGC 6946 (tarsenet "Zl liooih . The correlation between core 
radius and wa velength is also in the same sense as that found 
for M82-F by iMcCradv et al.1 (l2005h and may be an indication 
that mass segregation (primordial or dynamical) is present in 
NGC 1569-B. However, this needs to be confirmed by a more 
detailed analysis, including a modelling of how mass segregation 
would translate into observables such as core- and half-light radius. 

As a final note, we find it fascinating to contemplate that only 
a couple of decades ago, it was still uncertain whether NGC 1569- 
A and NGC 1569-B were in fact star clusters in NGC 1569 or mere 
foreground stars. Today, the Advanced Camera for Surveys on HST 
has made it possible to not only settle this question definitively, but 
to study the individual stars that make up these clusters. 
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